Longitudinal and transverse dissipation in a simple model for the vortex lattice with 

screening 
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Transport properties of the vortex lattice in high temperature superconductors are studied us- 
ing numerical simulations in the case in which the non-local interactions between vortex lines are 
dismissed. The results obtained for the longitudinal and transverse resistivities in the presence of 
quenched disorder are compared with the results of experimental measurements and other numer- 
ical simulations where the full interaction is considered. This work shows that the dependence on 
temperature of the resistivities is well described by the model without interactions, thus indicating 
that many of the transport characteristics of the vortex structure in real materials are mainly a 
consequence of the topological configuration of the vortex structure only. In addition, for highly 
anisotropic samples, a regime is obtained where longitudinal coherence is lost at temperatures where 
transverse coherence is still finite. I discuss the possibility of observing this regime in real samples. 



I. INTRODUCTION 

Physical properties of the vortex structure in high tem- 
perature superconductors are strongly dependent on the 
values of some external parameters, such as magnetic 
field, anisotropy, and quenched disorder .El Some of these 
properties -as for example the existence of a first order 
melting transition in clean samples-0 are originated in the 
interactions between vortices. Some others, however, are 
rather independent of the details of the interaction, and 
are related to the topological configuration of the vortex 
structure. It is interesting then to analyze in the simplest 
way the origin of those properties that depend only on 
the geometrical configuration of vortices, and not on their 
interaction. I concentrate in this paper on the behavior 
of the linear resistivity p (p = lim/_>o V/ 1) as a function 
of temperature -which has been widely used experimeru 
tally to unravel the properties of the vortex structure-, a 
both perpendicular and parallel to the applied magnetic 
field of a model that disregards non-local interactions be- 
tween vortices. The results of numerical simulations on 
this model compare well with experimental results ob- 
tained in YiBa2Cu3C>7 (YBCO, rather low anisotropy) 
and Bi 2 Sr 2 CaiCu208 (BSCCO, high anisotropy) sam- 
ples, as long as we consider zones of the phase diagram 
of the materials in which a first order transition is not 
observed. As stated above, the first order transition is 
generated by the interaction of vortices, and cannot be 
expected to occur in a model with only local interactions. 

I describe the model in the next section, and the results 
in section III. The case of very high anisotropies deserves 
special attention and the limit of two-dimensional sys- 
tems is discussed in section IV. The relevance of these 
results to real materials (in which vortices interact at fi- 
nite -usually large- distances) is discussed in section V. 
Finally, in section VI, I summarize and conclude. 



II. MODEL 

Vortices are modeled at different levels of detail when 
performing numerical simulations. A quite precise de- 
scription is the Ginzburg-Landau theory, formulated in 
terms of the superconducting order parameter W = 
1^1 exp (i#).el In this context, when an external magnetic 
field or thermal fluctuations are introduced, vortices ap- 
pear in the system as line-singularities around which 
§ 9 (r) dr = 2tt. A usual simplification which is appropri- 
ate in high temperature superconductors is to consider 
the modulus \^S\ of the order parameter as a constant, 
and keep only the phases 8 (r) as the dynamic variables. 
This leads to the study of the uniformly frustrated XY 
model,Q which has been extensively studied, both because 
of its intrinsic properties and because of its applications 
to superconducting systems. As first shown by Villain, □ 
the most important degrees of freedom of this model can 
be identified with the positions of the vortex in the sys- 
tem and an alternative description with these positions 
as the fundamental dynamical variables can be obtained. 
The original structure of the Ginzburg-Landau free en- 
ergy is reflected at this level in. ±he particular form of the 
interactions between vorticespQ which are cut off at dis- 
tances of the order of the some penetration lengths X a b, 
A c (the subindexes refer to the crystalline directions). 
These distances are in general large compared with the 
intervortex distance, and so the energy of the system has 
contributions coming from interaction between vortices 
that are far away from each other. I will consider here 
the case in which A is very small or, stated in another 
form, when the non-local terms of the interaction are 
dropped. In this case, the energy of a given configura- 
tion is simply proportional to the total length of vortices 
in the system. The investigation of the properties of the 
model in this case is important since, as we shall see, 
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it gives insight on the behavior of the system with the 
full interaction, and allows to understand that the origin 
of many properties of the vortex lattice comes from the 
topological structure of vortex lines, and not from the 
exact nature of the interactions. 

I will consider vortex segments lying on the bonds of a 
cubic mesh with periodic boundary conditions. Formally, 
the Hamiltonian of the system is 



(1) 



i,fj 



where are integer variables defined on the nodes i of 
the lattice, with direction p (// = a, b, c). The prime in 
the sum symbol indicates that only those configurations 
with zero divergence of the vector field n have to be con- 
sidered. The constants e,^ are the energies of vortex seg- 
ments at the positions i, p. We allow for the existence of 
anisotropy, defining a parameter r\ as 77 = (e^ c ) / (£i,ab) 1 
with (...) indicating averaged values throughout the lat- 
tice. Disorder is introduced by allowing the values of 
Ei tfi to be different in different points of the sample, fluc- 
tuating around the mean value. A disorder parameter 
(the same for the three spatial directions) is defined as 



D 



with e™ x and £™ n 



max pinirA / ( pinax _i_ _min 

being the maximum and minimum value of the energy 
of a vortex segment. The distribution between £™ n and 
e™^ x is taken flat. 

As the initial configuration of the system, a set of 
straight vortex lines directed along the c direction and 
uniformly distributed on the ab plane is considered. The 
number of vortices divided by the number of elementary 
plaquettes of the system perpendicular to the c direction 
defines the dimensionless magnetic field H. The Monte- 
carlo process for updating the configuration consists in 
sequentially proposing the creation of elemental squared 
loops in all plaquettes of the lattice and with the three 
possible directions. The acceptance of the new configura- 
tion is carried out using a standard Metropolis algorithm. 
The initial configuration and the Montecarlo procedure 
guarantee that at any moment the vortex configuration 
has zero divergence. In order to calculate resistivities, 
both parallel and perpendicular to the applied field (p c 
and p a b, respectively), we have to include a small ex- 
ternal current /. This is done by adding a term to the 
Hamiltonian (Q) that changes the energy of loops oriented 
perpendicularly to the current. One orientation increase 
its energy by +/, and the other one decreases it in the 
same quantity. The value of the external current / is 
chosen in such a way that the disbalance between right - 
and left-handed loops is never higher that 1/100 of the 
energy of the loop. In this regime of low currents the 
response is linear in the applied current, and resistivities 
do not depend on the exact value of I. The numerical 
results for different values of the anisotropy follow. 



III. RESULTS FOR THREE DIMENSIONAL 
SAMPLES 

I will describe the results of the numerical simula- 
tions performed in systems with progressively higher 
anisotropies. Three different regions are distinguishable. 
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FIG. 1. (a) Circles: resistivities along the ab and c direc- 
tions using the Hamiltonian (Q) for an isotropic lattice of size 
40 x 40 x 8 with D = and H — 1/4. The dotted line is a 
fitting to a thermally activated movement of single vortices of 
the form p a t ~ ^ exp (-2£o/teT). Stars: percolation proba- 
bility as defined in text, note that this magnitude is (different 
from) zero only if p c is (different from) zero, (b) The resistiv- 
ity values plotted to show thermally activated behaviors. 

For low anisotropies, typical results for p c and p a b as a 
function of temperature are shown in Fig. [I] for a sample 
of size L a x Lb x L c = 40 x 40 x 8 with D = (all £j iM 
equal to £ ) and H = 1/4. As we see the ab dissipation 
has a thermally activated behavior (dotted line fitting) as 
long as the c axis dissipation is zero (T < Tp). The ac- 
tivation energy of this process is 2£o, and corresponds to 
the energy necessary to create a double kink on a straight 
vortex, which is the first step in the process of movement. 
So in this region the dynamics of the system is that of 
a set of individual vortices thermally wandering into the 
sample. At a certain temperature Tp, the c axis dissipa- 
tion p c becomes finite. The origin of this dissipation is 
related to a percolation transition of vortices occurring 
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in the system, as studied before in. the three dimensional 
Josephson junction array model. u For T < Tp the vor- 
tices are practically independent, and a c axis directed 
current exerts no net force on them, thus generating no 
dissipation (in the linear regime). For T > Tp the vortex 
lines are so heavily interconnected that vortex paths run- 
ning (on average) in the ab direction appear. An external 
current applied along the c direction, being perpendicular 
to these paths, exert a net force on them and generates 
dissipation. In Fig. [l](a) we also see the value of the 
percolation probability P, defined as the fraction of time 
in which at least one percolation path is found in the 
system. It is clearly observed that the c axis dissipation 
is zero if P is zero. The transition in the P (T) curiae 
between and 1 becomes sharp in the limit L a b — > oon 
An interesting effect of the percolation transition on 
the values of p a b is observed. When p c starts to be dif- 
ferent from zero, the values of p a b deviate from the pre- 
diction of an activated behavior of individual vortices, 
and become smaller. This is an indication that for tem- 
peratures greater that Tp the percolation of the vortex 
structure and the existence of many thermal excitations 
interferes with the thermally activated behavior of sin- 
gle vortices. In fact, for T > Tp the concept of an iso- 
lated vortex in the system looses its sense, because we 
have a strongly entangled configuration of vortex lines. 
This effect of the percolation transition on the p a b {T) 
curve causes this to develop a typical shoulder that has 
been experimentally observed in measurements on YBCO 
samples, BE3 although its origin had not been clearly es- 
tablished. 

In the case D = 0, when increasing anisotropy, the ab 
plane activation energy goes to zero as 1 / 77 because this 
is mainly determined by the energy necessary to create 
a double kink, which is precisely 2/77. Also the c axis 
transition is governed by an energy scale of the order of 
~ I/77, because this temperature is mainly determined 
by the typical energy of the interlayer excitations (which 
goes as ~ 1 / 77). On heating, and from a practical point 
of view, the c axis transition occurs at a temperature at 
which p a b is clearly different from zero for all anisotropics. 

The physics is richer in the case of samples with defects 
(D ^ 0). In this case, when anisotropy is increased, the 
activation energy for the ab dissipation tends to a value of 
the order of D 1 / 2 , because the energy of a vortex segment 
piercing the ab plane is not constant in this case but has 
a dispersion of the order of D 1 / 2 , and this is the typical 
energy barrier that has to be overcome when a vortex seg- 
ment wanders within the ab planes. Anisotropy decreases 
the percolation temperature Tp in a factor ~ 1 /rj, but has 
minor effect on the thermal activation in the planes as 
long as D 7^ 0. Thus for high enough anisotropics, the c 
axis dissipation is expected to occur even at lower tem- 
peratures than the ab plane dissipation. However, when 
this range is approached a particular transition occurs in 
the system. The low temperature configuration of the 
vortex structure passes from a disentangled, and rather 
ordered configuration of vortex lines for low anisotropics, 



to an entangled configuration of vortex lines, i.e., we can 
say that the percolation temperature of the system drops 
abruptly to zero. The origin of this particular transition 
is the following. If anisotropy is low the system will pre- 
fer to remain ordered, with vortices almost straight in 
order to minimize their line energy. If anisotropy is in- 
creased, entangled configurations (in which vortices use 
the strongest pinning sites on the ab planes) diminish 
their energy and become locally stable, and from some 
value of the anisotropy, one entangled configuration be- 
comes globally stable. This is the critical anisotropy 771. 
This transition-is related to the Bragg glass to vortex 
glass transition^ proposed to occur in very anisotropic 
samples when increasing the magnetic field, because - 
as it has been discussed elsewhere-E3 an increase in the 
magnetic field is equivalent to an increase in the effective 
anisotropy and disorder of the system. 
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FIG. 2. Typical results for p c and p a b as a function of tem- 
perature for two different values of anisotropy, in a sample of 
size 40 x 40 x 8, with D = 0.5 and H = 1/4. 

Numerical simulations indicate that in the 7/ > 771 
regime both ab plane and c axis dissipation have a ther- 
mally activated behavior as shown in Fig. |^. The re- 
sults for activation energies of p c and p a b as a function of 
anisotropy for a system of 40 x 40 x 8, with D = 0.5 are 
shown in figure ^ (results for other values of H are sim- 
ilar, with only a rescaling of the anisotropy axis). From 
a numerical (or experimental) point of view, it must be 
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kept in mind that a thermally activated behavior with 
a given activation energy can be checked only for val- 
ues of UbT greater than some fraction (which depends 
on sensibility) of the activation energy. In particular, no 
statements can be made about the possible existence of 
a true critical temperature much lower than that corre- 
sponding to the activation energy. In Fig. [| and for 
f] < Vi ~ 4 the c axis activation energy is not defined, 
because the transition is a percolation process and not a 
thermally activated process, as discussed before. 
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FIG. 3. Activation energies for the dissipation along c and 
ab direction as a function of anisotropy (size 40 x 40 x 8, 
D — 0.5, H = 1/4). For n < r)i the c axis dissipation is not 
thermally activated but has a critical behavior instead. For 
n > rj2 the c axis activation energy decays as 1 jr\. 

The most direct and important conclusion from this 
figure is that there are two sub-regions in the range 
77 > r\\. For r\\ < r\ < r\2 activation energies for trans- 
port parallel and perpendicular to the field are rather 
similar, whereas for r\ > 772 the c axis activation energy 
is lower than the corresponding to the ab plane. In the 
whole range 77 > 771 , ab plane activation energy is rather 
anisotropy independent, indicating that the ab plane dis- 
sipation is governed by the thermal activation of ver- 
tical segments of vortex lines pinned to the ab planes. 
For 77 < 772 both parallel and perpendicular dissipation 
have the same activation energy indicating that they both 
originate in the same physical process S3 In fact, ab plane 
dissipation is caused by the thermal depinning of vortex 
lines that cross the sample in the c direction. In turn, 
the c axis dissipation is caused by the thermal depinning 
of vortex lines that cross the sample in the ab direction 
(which exist because the vortex configuration is heavily 
entangled). However, as the formation of these paths is 
mediated by the externally generated vortices they are 
also pinned to the ab planes, and the activation energy 
for both processes is the same. Despite this, a global fac- 
tor in the resistivity relation p c / p a b depending upon the 
anisotropy and the geometrical configuration of vortex 



lines (especially on the relation between number of paths 
running along c axis and ab direction) is expected. 

For 77 > 772 the decrease of the activation energy for c 
axis dissipation as 1 / 77 indicates that a new dissipation 
mechanism that depends only on the excitation of hor- 
izontal loops between planes is taking place. Being the 
vertical segments of vortex lines still frozen in the range 
of temperatures in which the c axis dissipation starts 
to be appreciable, this mechanism has to be related to 
processes occurring between consecutive ab planes, i.e., 
in the zone 77 > 772 a complete decoupling of the planes 
takes place. 

IV. PERPENDICULAR DISSIPATION IN TWO 
DIMENSIONAL SYSTEMS 

As I mentioned in the previous section, the activa- 
tion energy of the c axis dissipation at high anisotropics 
(77 > 772) goes to zero as I/77, indicating that processes 
involving only the horizontal segments between consecu- 
tive planes arc important. In order to understand clearly 
this kind of processes I will analyze the dissipation in the 
case of a unique horizontal plane. 

Consider the model studied previously (Eq. [I]) but 
now on a two dimensional geometry It is useful to point 
out that this model (in the case without disorder) has 
been naped-Jn other context the roughening model of 
surfaces. EJ~ll3 The mapping is made between vortex seg- 
ments in the original model and height differences of a 
growing surface in the roughening model, the divergence 
free condition in the vortex model being the key prop- 
erty that makes possible the mapping. In the language 
of the roughening problem, the surface is smooth at large 
scales in the low temperature phase, and is rough in the 
high temperature phase. The long extending terraces 
that make the surface rough are no more than the in- 
finite length vortex paths of the vortex model, and the 
growing rate of the surface maps onto the perpendicu- 
lar resistivity of the vortex model. In the homogeneous 
case (without disorder), the roughening model is known 
to have an inverted Kosterlitz-Thouless transition at a 
temperature Tkt — 0.36e where e is the energy of an ele- 
mental growing step £&n elemental vortex loop in the vor- 
tex model, e = 4eo)-li3 This temperature can be viewed 
again as a percolation temperature of vortex lines in the 
vortex system: below the critical temperature Tkt there 
is no infinite length paths, whereas for T > Tkt these 
paths exist. In Fig. |] the percolation probability for 
systems of different sizes is shown, and the transition is 
clearly distinguishable. As in the three dimensional case 
the transition in the variable P (T) becomes sharp in the 
limit L_b — > 00. From known results on the roughening 
modelEj we can directly conclude that for an infinite sys- 
tem the perpendicular resistivity (p c ) jumps from zero to 
a finite value at the temperature Tkt- 
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tem without disorder for different system sizes. A phase tran- 
sition in the limit of very large systems is clearly observed. 

The disorder in the two-dimensional case may have 
two origins: one is the disorder introduced directly in 
the Hamiltonian, through the dependence of Ei ifi on co- 
ordinates (Eq. (|])). The other comes from the existence 
in the 2D system of some quenched horizontal segments 
that come from the horizontal parts of the 3D vortices. 
This segments within the 2D system do not satisfy the 
condition Vn = 0, so its inclusion in the dynamics makes 
a non trivial contribution. We will concentrate on the 
second type of disorder for two reasons: numerical sim- 
ulations show that the qualitative changes produced by 
the two types of disorder are similar, and secondly be- 
cause this type of disorder is dynamical, and may change 
when changing anisotropy. 

So the problem may be stated as that of a Hamiltonian 



H = 



(2) 



where now bi ^ represent the horizontal segments induced 
by the 3D vortices, and n = a,b. Note that the variables 
rii^ still satisfy the condition Vn = 0, and the bare en- 
ergy of the segments lying on the links has been taken 
to be £0 in all sites. It is interesting to consider the be- 
havior of the system when the number of horizontal seg- 
ments is being increased. We can characterize this value 
as the fraction of links in which bi^ — ±1, (which we 
denote by D). A necessary condition for a finite dissipa- 
tion (perpendicular to the plane) is still that the vector 
field rii tfl generates a path running all across the (two- 
dimensional) system. The probability of existence of such 
paths is shown for different sample sizes and different 
values of the disorder in Fig. ||(a). For D < D cr ~ 0.3 
a well defined transition when the system size increases 
is detected. For finite disorder the transition tempera- 
ture Tp decreases with respect to the zero disorder value 
1.47. If disorder is too high, however, 



(D=0) _ rp 

there is no intersection of lines corresponding to differ- 



ent sizes, indicating that the percolation transition moves 
down to aero temperature. The behavior of the perco- 
lation temperature as a function of D is plotted in Fig. 
|(b). 
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FIG. 5. (a) Percolation probability in a two-dimensional 
system for different values of the disorder for system sizes 
20 x 20, 40 x 40, and 80 x 80 (from bottom to top, at the right 
of the curves). The results are the averaged values over 3 con- 
figurations of disorder for D = 0.125, 0.25, and over 8 config- 
urations for D = 0.375, 0.5. For D < D cr ~ 0.3 a percolation 
phase transition at finite temperature is clearly detectable. 
For higher values of disorder percolation temperature moves 
down to T — 0. (b) Temperature of the percolation transition 
as a function of disorder as obtained from the results in panel 
(a). Points are the results of numerical simulations. Dotted 
line is a guide to the eye only. 

Perpendicular resistivity simulations in this model 
(Fig. |^) give values that cleanly go to zero in the case in 
which Tp is finite (D < 0.3)lia, and in the case when this 
variable is finite at any temperature give results that can 
be well fitted by a thermal activation expression. How- 
ever, the possibility that the system has a real critical 
temperature at a temperature where the resistivity be- 
comes undetectable in the simulations cannot be ruled 
out (see the note on the mapping to the roughening prob- 
lem with disorder in the next paragraph). 
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FIG. 6. Perpendicular resistivity of a two dimensional sys- 
tem (80 x 80) as a function of temperature for different 
values of the concentration of quenched segments D. For 
D < D cr ~ 0.3 the approximate temperature at which tem- 
perature vanishes is indicated by arrows. For D > D cr the 
curves show activated low temperature tails. 

The behavior of the percolation probability and resis- 
tivity for different values of D of model (0) can be qual- 
itatively explained using an argument of the type of the 
celebrated original Kosterlitz-Thoulcss argument for the 
phase transition in the two-dimensional XY model. It 
goes as follows. The energy of a path running across a 
system of size ~ L is the sum of ~ L variables. According 
to (g), each of these variables may have the value e , —£o, 
or 3eo depending if for that particular size the variable 
b is 0, ±1. The energy of the path becomes a Gaussian 
variable with mean value ~ L and dispersion y 'LDeq. 
The total number of paths of length ~ L running across 
the system is exponential in L, i.e., about ~ e L . Taking 
into account these estimations, and considering the calcu- 
lated number of paths with different energies as a density 
of states for non-interacting "particles" , it is possible to 
make the statistical mechanics of the system and look for 
the thermodynamic free energy of the equilibrium config- 
uration. I only show the results, and not the details of 
the calculations, which are straightforward. Two clearly 
different regions appear (see Fig. ^): there is a zone (la- 
belled I) of high disorder or high temperature in which 
some percolation paths have negative free energy, and 
so their existence is thermodynamically stable. In the 
other zone (labelled II) with low temperature and disor- 
der, all percolation paths have positive free energy and 
thus the configuration without any path is the stable one. 
In zone II the system does not contain any percolation 
paths, and its resistivity is zero. When passing to zone I 
increasing the temperature, many percolation paths be- 
come immediately accessible at the transition, and per- 
pendicular resistivity jumps to a finite value. In the case 
D > D cr there are available percolation paths down to 
zero temperature, and its likely that resistivity remains 
finite (with a thermally activated behavior) up to zero 



temperature, although a critical behavior of the resistiv- 
ity at some finite temperature cannot be. completely ruled 
out from these qualitative arguments.lI3 These estimates 
compare well with the numerical results on Hamiltonian 
(§) (Figures @ and (§). 
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FIG. 7. Qualitative phase diagram from the estimation of 
free energy of percolation paths, as discussed in the text. Zone 
I is a region with finite linear dissipation. Zone II has zero 
linear dissipation. The labels la and lb and the dashed region 
are linked to the discussion presented in [18]. 

The mechanism discussed in this section for the per- 
pendicular dissipation of a two-dimensional system is the 
responsible of the c axis dissipation in the three dimen- 
sional case for 77 > 772 discussed in the previous sec- 
tion. As it was indicated before this mechanism relies 
on the movement of vortex segments between consecu- 
tive ab planes, contrary to the dissipation occurring for 
77 < rj2 which involves paths that occupy a finite thickness 
along the c direction. 

V. RELEVANCE TO REAL MATERIALS. 
EFFECT OF NON LOCAL INTERACTIONS 

In spite of its simplicity, the previously studied model 
reproduces many characteristics of real transport mea- 
surements in high-T c superconductors, indicating that 
some of these properties are rather independent of the 
exact form of the interaction, and they are only conse- 
quences of geometrical properties. In particular, the case 
i] < r/i (Fig. |3|) gives results "wiji^h. are comparable to ex- 
perimental results in YBCOJl9'c2l In this case the c axis 
resistivity has been observed to appear only at higher 
temperatures than the ab dissipation, and a shoulder in 
Pab near Tp similar to that of Fig.|l| has been experimen- 
tally observed. Numerical simulations using the complete 
vortex-vcptex interaction give also the same qualitative 
behavior OI The results of the previous sections for the 
activation energy in the range r\\ < 77 < 772, namely the 
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coincidence of activation energies for p ab andpOr, has been 
experimentally observed in BSCCO samplestj and theo- 
retically interpreted as an indication, of the same origin 
for both dissipations in this regime£3 

The only one region that gives a qualitatively new re- 
sult is the range 77 > 772 , where a finite value for p c was 
obtained even at temperatures at which p a b is still negli- 
gible. We have to determine if this region can occur when 
the complete vortex-vortex interaction is taken into ac- 
count or if it is a consequence of the local form of the 
interaction used in the model. The complete Hamilto- 
nian for interacting vortex segments can be written 
in the generic form 

H = Y1 G^{Y t -Y j )n^{Y i )n^(Y J ) . (3) 

where G M (r) is the non-local interaction function be- 
tween segments (in previous sections I used this type of 
Hamiltonian with a local function G). The function G ao 
(G c ) is globally proportional to an energy scale that we 
can identify with our e ab (e c ) of equation (||), and the 
naive conclusion would be that in the highly anisotropic 
case (e a b *C s c ) the system of interplane vortices has 
a transition temperature T% of the order ksT\ ~ e a b, 
whereas the system of vortex-antivortex pairs within the 
planes has a transition temperature T2 of the order ksT^ 
~ s c 3> ksTi. However, the structure of the function G 
for segments parallel to the layers 

G ab (k) = 

_ Sab 

~ [4-2 cos (k x ) - 2 cos (k y )} + [2 - 2 cos (k z )] + ^ 

(4) 

(where \J A/e a h = A c is the magnetic penetration 
depth along the c direction in units of the lattice spac- 
ing) shows that the range of the interaction of paral- 
lel to plane segments increases with anisotropy as ~ 
(min [e c , A] /eq},) 1 ^ 2 , and this cannot be neglected EH 
In fact, in the language of the renormalization group 
and in the case of very large A c , interacting horizon- 
tal segments renormalize to non-interacting segments but 

with an energy ~ e a b {(ec/sab) 1 ^ 2 j = £c, so the tran- 
sition temperature T\ is of the order of ksTi ~ e c . A 
more detailed calculation! gives Tj ~ Am min [e c , 2A], 
T2 ~ ? min [e c , 2A] . To have finite p c at temperatures 
where p ab is still zero is necessary at least that T\ be 
lower than T 2 , and this is not the case. 

These estimates seem to rule out the practical possi- 
bility of a zone like that for r\ > r\i in figure |^. However, 
the analysis made corresponds to the case of zero exter- 
nal field and disorder, and we know that the inclusion 
of them decreases the transition temperature of the in- 
terplane system of loops. So the question if in a finite 
external field, and for high anisotropics we can have a c 



axis transition temperature lower than the corresponding 
to the ab plane has still to be answered. The point cannot 
be fully discussed in all situations, but there is a limiting 
case in which it can be addressed. In fact, the case with 
a divergent anisotropy, in which the horizontal segments 
interact through a really long range potential, can sur- 
prisingly be discussed much easily that the case in which 
the interaction is local. Let us consider a finite value of 
e a b, and take e c , A — > 00, (so T2 — > 00 also), and analyze 
the problem of the system of horizontal vortex segments 
between two consecutive planes. To be able to work with 
some indeterminations that will appear, the system size 
we will consider the system size will be considered finite 
with a value L x L, and the limit L — > 00 will be taken 
at the end. In this limiting case the function G a b (k) in 
terms of the two dimensional vector k = (k x ,k y ) , reduces 
to Gab (k) = e a b/ [4 — 2 cos (k x ) — 2 cos (k y )] . Introduc- 
ing (integer) loop variables 1(y) through the substitution 
d^l (r) = (r) , and integrating by parts twice in the 
Hamiltonian, an equivalent model is obtained: 

H = e ab Y,(l(r)-T) 2 + lJ2 l (r), (5) 

r r 

where I = ^2 r l (r) /L 2 , is the mean value of the loops 
variables over the L 2 plaquettes of the system. The time 
derivative of / is proportional to the voltage generated 
by the external current /. This expression is much easier 
to handle than the original formulation, and it is a sort 
of mean field Hamiltonian in which variables I interact 
only through the I term. I stress however that no ap- 
proximations other than the ones mentioned were done 
on passing from (^) to (||) . When L is large, the partition 
function corresponding to (^|) can be calculated using a 
Lagrange multiplier for I. The energy as a function of I 
becomes a periodic function of I (except by the current 
term) , and in the limit T s ab it reads 

E = L 2 {- r ^- e ~^ T/eal cos ( 27r + ■ ( 6 ) 

When L — > 00, and for any value of temperature there 
is a critical value of /, namely I cr — T / £ a& _ This 

critical current is nonzero for any temperature and thus 
indicates that the critical temperature of the system is 
infinite. This is consistent with the exact results since we 
took e c — > 00. Notice, however, that the extremely small 
value of this critical current when T is much larger that 
e a b turns it difficult to verify this results in numerical 
simulations. 

When quenched disorder due to horizontal segments 
of externally generated vortices are included, the previ- 
ous picture changes in the following way. The poten- 
tial due to these quenched segments is also long ranged, 
and on the / variables they produce a potential that 
goes as ~ s a b/r. Summing up on each site the con- 
tributions from a random distribution of disorder on a 
sample of size L gives a potential of typical amplitude 
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V ~ e nh D 



1/2 



£ a bD\j2ix In L, where D 

is the concentration of quenched segments. The Hamil- 
tonian now reads 



r r 

+ e ab DV2nhiLj2xrl (r) . (7) 

r 

The Xr are random variables with typical value 1 and 
satisfying ^ r \r = Oo In the same way as before the 
energy of the system may be written in the form 

E D = L 2 (-Tl e -^ T ^»e-^ D2 111 L ' 2 cos (2ttZ) + jf) , 

(8) 



E D = -^—e-* T ^L 2 -« D I 2 cos (2ttZ) + L 2 /F. (9) 

Sab 

When L goes to infinity the first term will give an in- 
finite activation energy -and thus a zero resistivity- if 
D < 2/tt 3 / 2 . On the other hand if D > 2/tt 3 / 2 acti- 
vation energy goes to zero and the resistivity is finite. 
D cr = 2/tt 3 / 2 is the critical value of the disorder in the 
system. As we see, in this case the c axis transition occurs 
at zero temperature, whereas the ab transition tempera- 
ture is ~ e c — > oo. 

Thus at least in the case of infinite ranged interactions 
between vortex loops, it can be analytically shown that 
the c axis transition occurs at lower temperature than the 
ab transition, if disorder is greater than a critical value. 
This suggests that even in cases with finite e c , disorder 
can make longitudinal resistivity to appear at lower tem- 
peratures than transversal dissipation. The conclusion is 
that with the inclusion of the full interactions between 
vortices and for finite values of e a b and e c (e c 3> e a b), a 
phase diagram as the one sketched in Fig. |?] still holds, 
but with the zero disorder temperature transition being 
proportional to e c rather to e a b- This indicates that a 
regime where p c ^ and p a b = may be experimentally 
accessible in highly anisotropic and disordered samples. 



VI. SUMMARY AND CONCLUSION 

In this paper I have presented results from numerical 
simulations of vortex lattices in the presence of exter- 
nal magnetic field and quenched disorder, in which the 
interactions between vortices are neglected. In spite of 
the simplicity of the model the results reproduce qualita- 
tively well many characteristics of the p c (T) and p a b (T) 
functions observed in experiments, both for YBCO sam- 
ples, in which p c becomes different from zero when the 
value of pab is already appreciable, and for BSCCO where 



experiments indicate that p c (T) and p a b (T) have ther- 
mally activated behaviors with the same activation en- 
ergy. In addition, a new regime in which activation en- 
ergy for p c (T) is lower than that for p a b (T) was found 
in the simulations, corresponding to the case where co- 
herence length along c direction reduces to the distance 
between consecutive superconducting planes. Although 
in principle this regime may be an artifact of having dis- 
regarded interactions, I have shown it may occur in real 
samples in the presence of disorder, external magnetic 
field, and for high enough anisotropy. 
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